Quantum Hydrodynamics of Fractional Hall Effect: 
Quantum Kirchhoff Equations 
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We argue that flows of the quantum electronic liquid in the Fractional Quantum Hall state are 
comprehensively described by the hydrodynamics of vortices in the quantum incompressible rotating 
liquid. We obtain the quantum hydrodynamics of vortex flow by quantizing Kirchhoff equations for 
vortex dynamics. We demonstrate that quantized Kirchhoff equations capture all major features of 
FQH states including subtle effects of Lorentz shear force, magneto-roton spectrum, Hall current 
in a non-uniform electromagnetic field, thus providing a powerful framework to study FQHE and 
superfluids. 
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1. Introduction. In Fractional Quantum Hall (FQH) 
regime electrons form a peculiar quantum liquid. Charac- 
teristic features of the liquid are: flows are incompress- 
ible [lj], almost dissipation- free 0, Q, the Hall conduc- 
tance is quantized Q, vortices are elementary excitations, 
vortices carry fractionally quantized negative electronic 
chargesjpl, they are separated from the ground state by 
a gap [J, |j| . More subtle features are the universal an- 
harmonic term of the structure factor and magneto-roton 
minimum in excitation spectrum [3] (see ([M]) ), quantized 
double layers of a density and shear at boundaries and 
vortices the Lorentz shear force @. 

Such liquid does not posses linear gapless waves except 
edge modes propagating along the boundary. Only avail- 
able bulk flow is a non-linear flow of vorticity. Since in 
FQH regime vorticity is linked to the electronic charge, 
flows of vorticity are related to a charge flow |25j . 

A natural approach to study flows of FQH states is 
hydrodynamics. Not all quantum liquids are subject of 
quantum hydrodynamics. Quantum hydrodynamics is 
based on a fundamentally restrictive assumption that all 
quantum states are fully characterized by the density and 
velocity. Apart from of superfluids (and superconduc- 
tors) and Luttinger liquids, FQH liquid is, yet, another 
important case. 

Hydrodynamics of quantum fluids goes back to Landau 
Q and Feynman Q, but remains in its infancy. A quest 
for the hydrodynamics of FQH liquid has been originated 
by a seminal paper Q , further discussed in [9H12| , it is a 
focus of a renewed interest. 

As always in hydrodynamics only a few basic princi- 
pals, symmetries, and a few phenomenological parame- 
ters are sufficient to formulate fundamental equations. 
To this aim an underlying microscopic Hamiltonian de- 
scribing emergent FQH states under a strong Coulomb 
interaction is in fact not necessary. 

In the letter we attempt to formulate these principals 
and to develop hydrodynamics of FQH bulk states in a 
close analog of Feynman theory of superfluid helium |8[ 
and magneto-roton theory of collective excitations by [J] . 

We argue that flows of FQH liquid are equivalent to 
flows of vortices in the quantum incompressible rotating 



Euler liquid. Based on this conjecture we obtain the ma- 
jor features od FQHE including subtle effects of Lorentz 
shear force [f|and magneto-roton spectrum H missed by 
the previous hydrodynamics approaches [9Hl2j]. 

Hydrodynamics of vortex flows (classical and quantum 
alike) is an interesting subject by its own. Apart from 
FQHE it is also relevant to the theory of superfluids and 
classical hydrodynamics. 

Quantization of hydrodynamics is a subtle matter. 
It is best achieved through quantization of Kirchhoff 
equations [13]. Classical Kirchhoff equations describe 
a motion of vortices in a 2D incompressible isentropic 
fluid. We show that, quite remarkable, quantum Kirch- 
hoff equations capture all known features of FQH liquid. 
They can be used as a platform for studying FQH Effect. 

In this paper we consider only the simplest Laughlin's 
cases where fraction v is an inverse of an odd integer, 
say 1/3. In these cases electronic liquids do not posses 
additional symmetries. Extension to other FQH states, 
not discussed here, is possible and interesting. 

We start by reminding classical Kirchhoff equations for 
rotating incompressible Euler flows (see e.g..[14|). 

2. Classical Kirchhoff equations. 2D incompressible isen- 
tropic flow (that is a flow where gradient of density are or- 
thogonal to gradients of pressure, or zero) is fully charac- 
terized by its vorticity. Vorticity obeys a single (pscudo) 
scalar equation, which in the case of inviscid fluid has a 
simple geometrical meaning: the material derivative of 
the vorticity vanishes 
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u- V \w = 0. 



(1) 



In other words vorticity zu = V x u is transported along 
divergent-free velocity field u 



V-u = 0. 



(2) 



Helmholtz, and later Kirchhoff realized that there is a 
class of solutions of the vorticity equation (flj which 
consists of a finite number of point-like vortices mov- 
ing by Magnus forces. In this case the complex velocity 
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u = u x — iu y is a meromorphic function 

N r 

U ( Z ,t) = -mz + iJ2— 4777, (3) 
i=l * l ' 

where O is an angular velocity if the fluid is rotated, N 
is a number of vortices, Ti and £j(t) are circulations and 
positions of vortices. We denote z = x+iy, d= \{V x — iV y ) 
and use the roman script for complex vectors a = a x —ia y . 

A substitution this " pole Ansatz" into the Euler equa- 
tion ([]} yields that the number of vortices N and the 
circulations T t do not change in time, but moving posi- 
tion of vortices Zi(t) obey the Kirchhoff equations. 

If the rotation of fluid is very strong, vortices prefer 
to be of the same sign opposite to the direction of the 
rotation. Bearing in mind the quantum case we assume 
that vortices have the same (minimal) circulation 1^ = T. 
Then Kirchhoff equations obtained by localizing the Eu- 
ler equation to vortex cores at z = Z{ expresses velocities 
v, of vortex cores through their positions 

N r 

* s *=- m+ %itirwY (4) 

The Kirchhoff equations replace the non-linear PDE (JTJ 
by the dynamical system, reflecting integrability of the 
Euler flow. They can be used for different aims. Equa- 
tions describe chaotic motions of a finite number of vor- 
tices if N > 3. Alternatively they can be used to approx- 
imate virtually any flow if N — > oo. 

3. Chiral flow. The flow relevant for FQHE is the chiral 
flow, where a large number of background vortices uni- 
formly distributed with the mean density p — Q,/(ttT) 
largely compensates rotation. This is a very special flow 
in fluid mechanics. There we distinguish two types of mo- 
tion: fast motion of the fluid around vortex cores, and 
a slow motion of vortices. In this respect vortices them- 
selves must be considered as a (secondary) fluid. At the 
ground state of chiral flow vortices do not move, but the 
fluid does. Vortices carry an inertia m* - a phenomeno- 
logical parameter determined by the equation of state. It 
is not directly related to the inertia of fluid. 

In this paper we propose to model FQHE by a quan- 
tized slow flow of vortex fluid. 

4- Quantization of Kirchhoff equations. Quantization of 
the Kirchhoff equations consists of three steps: canonical 
quantization, a choice of the representation and the inner 
product. 

It is convenient to use a circulation q = m*r of mo- 
menta of vortices m*«j. The Poisson brackets followed 
from canonical Hamiltonian structure of the hydrody- 
namics. They are equal to the volume per particle per 
circulation {xi, yj}p.B. = l/farpq). We replace them by 
commutators 

{zi, zjjp.B. -> [Ei, Zj] = 2t 2 5ij. (5) 



Here 2£ 2 = v/(irp) has a dimension of area. Dimension- 
less number v = H/q is a semiclassical parameter. 

The next step is a choice of states. We assume that 
states are holomorphic polynomials of Then operators 
Zi are canonical momenta 

Zi = 2£ 2 d Zi . (6) 

The third step is to specify the inner product. We impose 
the chiral condition: operators Zi and Zi are assumed to 
be Hermitian conjugated 

chiral condition : Zj = zj. (7) 

This condition combined with representation ([6]) identi- 
fies a set of states with the Bargmann space [la]: the 
Hilbert space of analytic polynomials ip{z\, . . . , Zn) with 
the inner product (exp(— zz)dexp (zz) = z insures (|6I7[1 ) 

ty'ty) = J dp^ip, dp = Y[e~^ 2 ~d 2 z i . (8) 

i 

Consequently zj = —2£ 2 B Zi + z%. The normal ordering in 
the Bargmann space assumes that holomorphic operators 
stays to the left to anti-holomorphic operators. 

Eqs. (|4I6[) help to write velocity operators of vortices 

T- 1 v i = -i2vd Zi +i^ j —^—, (9) 

Eqs.((3]|ni) are quantum chiral Kirchhoff equations. They 
are readily generalized to a sphere, or a torus. 

5. Quantum Chiral Kirchhoff Equations and FQHE. We 
identify quantum chiral Kirchhoff equations with FQHE. 

First we comment that Bargmann space of analytic 
polynomials is another way to say that all states belong 
to the lowest Landau level [l|| (and also 0, H3|). Then 
we assume that a state where vortices are at rest is the 
ground state of the system ipQ. It nulls all velocity oper- 
ators Vj Arg [-0 O ] = (velocity operator acts on the phase 
of the w.f.). A common solution of the set of 1st order 
PDEs is the holomorphic part of the Laughlin w.f. 

i'o = l[(z l -z J f, hp = q. (10) 

i>j 

The w.f. becomes single valued if (3 is integer, antisym- 
metric if (3 is an odd-integer. 

In this interpretations vortices are identified with "par- 
ticles" entered into the Laughlin function. Hence a quasi- 
hole iph = rii( z o — Zityo [J is a hole in the uniform 
background of vortices - an anti-vortex. Thus we as- 
sign the electronic charge to vortices and identify an- 
gular velocity with the cyclotron frequency of vortices 
£1 — eB/(m*c). Then the entries of the Kirchhoff equa- 
tions are the magnetic length £ = yJKcJeB and the filling 
fraction v=phc/eB = h/q (T — n/(irp). 

The phenomenological parameter is the inertia of 
the vortex. It is naturally to assume that the energy 
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associated with the inertia A„ = Ml = h 2 /(£ 2 m*) is of 
the same order as the energy of a quasi-hole at rest, or 
equivalently, is a gap in the excitation spectrum. The 
latter is of the order of Coulomb energy e 2 /£. It is known 
experimentally ~ 10K [3j. The very existence of the 
FQH state requires that this energy scale must be less 
than the cyclotron frequency hu c ~ 25meV Ml, or 
that m* exceeds the bare electronic mass m* 3> m e . An 
assumption that the energy of the quasi-hole at rest is the 
gap and is related to the inertia of the flow is the major 
physical input led to the hydrodynamics of FQHE. 

A noticeable feature of the identification is that slowly 
varying external fields (potential well, gradients of tem- 
perature etc.) are coupled to slow moving vortices. 
For example a potential well U(r) adds to the energy 
X^i U{ri), where r< are coordinates of vortices, not fluid 
particles. It yields the Lorentz force —i[U, Zj\ = i2£ 2 d Zi U 
into the r.h.s. of (JSJ). In the flow which keeps the center 
mass at the origin, Kirchhoff equations yield £• Vj = 
i(irqp) ~~ 1 J2i 9 Zi U. This implies the fractionally quan- 
tized Hall conductance a = e 2 /(2nq) = u{e 2 /h). 

Finally we emphasize that the velocity (J5J) is different 
than velocity of individual electrons and coordinates of 
their guiding centers. 

In the rest of the paper we show that Kirchhoff equa- 
tions contain other, more subtle properties of FQHE. To 
this end we must develop the hydrodynamics of quantum 
vortex flow. To the best of our knowledge this has not 
been done even for the classical fluids. 

6. Velocity field of the flow of vortices. Eulerian hydro- 
dynamics of vortex computes the flux of vortices 

p = _^YJs(r-n),Vi} = -^{p,v}=m*^/pv^/p, (11) 

i 

where r,, are coordinates and velocities of vortices, {, } 
is the anti-commutator and p(r) — 8(r — r») is the 
density of vortices. This formula also defines the velocity 
of the vortex flow v(r). 

By construction operators P and P* are ladder opera- 
tors. They annihilate the ground state PipQ = V'oP^ = 0. 

In hydrodynamics the mass density of the fluid and 
mass flux of the flow are independent fields obeying the 
continuity and the Euler equations. Contrary, the ve- 
locity of the vortex flow v(r) (as well as the velocity of 
the fluid u(r)) are expressed through the density of vor- 
tices p(r). We must determine these relations. Before we 
proceed the comment is in order. 

In 2D Euler liquid the mass density of the fluid and 
mass density of vortices, both, obey the continuity equa- 
tions. Hence an initially imposed local condition, n(r) — 
m*p(r) is compatible with the dynamics at all time. This 
condition reduces the set of states being essentially equiv- 
alent to the chiral constraint (0. Thus the vortex flow 
can be formally considered as a special flow of 2D Euler 
liquid with a constituency relation between the flux and 
the density of the liquid [26j]. From this point of view 
J = \{n, u} = ^{p, u} is the mass flux of the special 



flow. It includes fast fluid motion around vortices and a 
slow motion of vortices. This is a useful but an auxiliary 
object. Instead we are after the flux of the vortex flow 
(fTTj) . which consists of only slow motion of vortices. 

The flux of the vortex flow subtly differs to the flux of 
the fluid. We compute them now. 

The velocity of the fluid is easy obtain. It is given 
by ([3]), where we represent the operator z by use 
pV-r- = S(r — rj)V rj , and pass to the continuum limit 
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5p 



(12) 



The flux J = ^{p,u} follows. Computing the flux of 
the vortex flow (fTl"]) we use: the 9-formula irS = d(-) and 
the identity 2 . = (£. -1-)*-^^-)*. 

With the help of dH) a simple computation yields 

i§3p = J + i§5p. (13) 



<1 



We obtain the important relation between flux of the 
vortex flow and the flux of the fluid. Using the notation 
a p\ = ([tv a v for 2-vectors we obtain 



P = J+iV'p, 
4 



v = u + — p V*p, 



(14) 



The shift (fT4j) holds in the classical and the quantum 
cases. It has far reaching consequences 27 1. 

The shift (fl4^) can be seen as a similarity transforma- 
tion. It preserves the volume. Hence the flow of vortices 
is incompressible V • v = like the fluid itself. 

The shift has a simple meaning. Velocity of the fluid 
u diverges at a core of an isolated vortex (as is in ([3])). 
However, velocities of vortices are finite. The shift re- 
moves that singularity. 

Further meaning of the shift is seen from monodromy 
of the wave function. Monodromy is the circulation of 
each particle in units of V. It is equal to the number 
of magnetic flux quantum piercing the system Nq. The 
circulation of a particle (a vortex) around the system of 
remaining N — 1 vortices is T(N — 1) = 2ir § v ■ dr. The 
monodromy is Nj, — /3(N — 1), i.e., the number of zeros 
of the w.f. with respect to each coordinate. On the other 
hand the circulation 2n <f> u ■ dr gives the total charge 
TN. The shift amounts for the difference. It simply 
means that a vortex does not interfere with itself. Eq. 
(fl"4]) can be seen as a local version of the global condition 
N<j, = fiN - 2s, where the shift hs =q/2 

We summarize the formulas for flux and velocity 



P = \{p^ P }-\pV*y 



<l 



|v*logp. 



(15) 
(16) 
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Potential <p obeys the Poisson equation Aip — -An(p— p). 
It is chosen such that flux vanishes at the ground state. 

7. Chiral constituency relation. The next step is to ex- 
press the flux of the vortices in terms of their density. 

Using the formula [p, dTT p ] — —\h~dp and the chiral re- 
lation 2l 2 d\. = ~2£ 2 d Zi + Zi, we write 

{p, dTT p } = -2dnjp + [p, dn p ] = -ih(d + -^—)p. 

Applying this formula to (|15l) we obtain the chiral con- 
stituency relation. We write it in two suggestive forms 

P = -pV*#, y = \[p-{\-v)\og P ], (17) 
P = [ -dT, T=\{d^-{\-v)d i t p. (18) 

TT Z Z 

The field VP has a meaning of the stream function, hence 
vorticity is 

uj = -A^=2nq[p-p+^-{^-u)Alogp]. (19) 

We observe that vorticity of the chiral flow differs from 
the density of vortices by the term oc A log p. The co- 
efficient comprises of the shift and quantum correction. 
The source of quantum corrections is the difference of the 
vorticity of the fluid m — 2irq[p — p— j^A log p] from the 
density of vortices. 

Integration of (|17l) gives the global version of the chiral 
condition. It connects the moment of inertia L = J{rx 
P)d 2 r and the gyration G = Jr 2 (p — p)d 2 r of the flow: 

£ 2 L = hG + Ne 2 h(/3 - 2) 

This form generalizes the exact sum rule of the Laughlin 
state at P = : Ei(0|N 2 |0) = Nl 2 (N — ((3 — 2)) . 

Eq. (fTT)) can be used to find density profiles for various 
coherent states_[28j. Consider e.g., a quasi-hole if>h = 
— Zi)if)o This state describes an anti-vortex (a 
hole in the sea of vortices). Its flux is P = ifi z f z ■ Eq 
(fTT)) expresses its density through the 2-points function. 
Eq.(??) computes the change of the gyration by a quasi- 
hole: N- 1 J r 2 Spd 2 r = -I 2 = This result often 
interpreted as a fractional charge of the quasi-hole - it 
occupies a fraction of volume per particle. Outside of the 
core, the quasi-hole as a source for vorticity in a classical 
equation [29[ 

-5(r-r ) =p[p-p+— (i-^Alogp]. 

8. Governing equation for the vortex flow. Since the den- 
sity of vortices determines its flux, the continuity equa- 
tion m*p + V • P = is the only governing equation of 
the chiral flow. Eq. (fTT|) provides the close form 

r 

p+ -Vip x Vp = 0. 



This is of course the standard equation for vorticity of 
incompressible liquid. Peculiar effects are burred on 
boundaries or in initial data. A moving edge must be 
a stream line (|17p . Since stream lines are non- linearly 
connected to level lines of the potential ip the dynam- 
ics differs from that of a free boundary of incompressible 
fluid. In [j| we showed that edge dynamics is governed by 
Benjamin-Ono equation. Similarly, initially given vortic- 
ity is non-linearly related to the density of vortices, hence 
its dynamics differs from Euler flows. 

9. Hamiltonian structure of vortex dynamics. Kirchhoff 
equations @ are integrable. Same equations follow from 
different Hamiltonians |3lj . One Hamiltonian comprises 
by kinetic energies of the fast fluid motion and slow drift 
of vortices He = I ^ p~ 1 3d 2 r. An equivalent form 
of the classical version of this Hamiltonian (but without 
incompressible and chiral conditions), has been proposed 
as a description of FQHE in 

Our main result is 

that flows of FQH states are described by another Hamil- 
tonian (used in Ql, see also [30j). containing only slow 
motion of vortices [H, HH 

H = — !— / P^P d 2 r. (20) 
2m* J 

Both yield the same Kirchhoff equations but have differ- 
ent energy. The difference is essential. The latter yields 
correct properties of the Laughlin state reflected by (|14[) , 
the magneto-roton minimum in the excitation spectrum 
of [J] and Lorentz shear force discussed below. The for- 
mer does not [13, HI] . 

With the help of the chiral relation (fTTj) we obtain 

H = [ p\V^\ 2 d\ (21) 

where the stream function is given by (|17|) . 

The symplectic structure can be expressed in two 
forms: in terms of velocity or flux of the fluid, or in terms 
of the vortex fluid. The first structure is the canonical 
Heisenberg current algebra followed from (fT2|) 

ml[u(r), u f (r')} = 4nq5(r - r'). (22) 

The second is an extended algebra obtained by the trans- 
formation (fT4"l) [36l W?\ . It is conveniently written in 
terms of the flux 

[P(r),P t (r')]= 2h[-Px\7 + qp(2np+jA)]6{r-r'). (23) 

The advantage of the first structure is its canonical form. 
The advantage of the second structure is that operators 
P , P^ are ladder operators annihilating the ground state. 

The Hamiltonian (|2"0)) and the algebra (|2"3"]l acting in 
the Bargmann space are the compact and comprehensive 
formulation of FQHE. 

10. Structure function and excitation spectrum. We il- 
lustrate the current algebra by computing the celebrated 
results of IJI for the structure factor and the excitation 
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spectrum. Let pk is a Fourier mode of a small den- 
sity modulation. In the linear approximation (|17[) con- 
nects it to the flux mode P fe = hj^il - ^-(k£) 2 )p k . 
Then project (f23|) onto the ground state and compute 
its r.h.s.. This gives the correlation of flux modes 
(P fe pt) = 2ph 2 l- 2 {l-^(M) 2 ). Comparing we obtain 

the Feynman-Bijl formula [8] for the single mode approxi- 

/p pt \ 

mation of the excitation spectrum A(fc) = 2m -I , and the 
structure factor s(k) — (pkP-k) 



A(k)- 



2 1,2 



K 2 k 



2m*s(k) ' 



S(k): 



2l-^(W) 2 ' 



(24) 



The excitation energy A(fc) = 2vr/3^-(l-^(M) 2 ) shows 
the magneto- roton minimum of [4|, the leading order of 
the structure factor also agrees with 

11. Stress tensor and Lorentz shear force. The shift (|14[) 
causes subtle hydrodynamic phenomena. One of them is 
the Lorentz shear force, a conservative force acting nor- 
mal to the shear flow. It enters to the conservative part of 
the viscous tensor, which we now compute. Initially be- 
ing introduced for the integer QHE in [19( and extended 
to FQHE in [2(J [Hj], the Lorentz shear force is in fact, 
the classical phenomena. It is a feature of vortex flow. 

We cast hydrodynamics equations in the form of the 
conservation law with the Lorentz force F = eE—-Bv * 



(25) 



and compute the momentum flux H^. 

We start from the first Hamiltonian structure writing 
conservation law for the fluid flux J: dt J M + V^II A1 , y =0. 
In this case the canonical property of the first structure 
(|22p determines the form of the momentum flux: IT^ = 
J l i( , m*p) 1 J» + PV^^v, while incompressibility condition 
determines the intrinsic pressure p (all operator products 
are normally ordered). 

We obtain the momentum flux of vortices by 
the similarity transformation (|14[) . Taking time deriva- 
tive of (|14[) and using the continuity equation, we write 



P—i 

M ~ 4 



V*(v ■V)p=V u (a flI/ -Il^). This gives the 



transformation of the momentum flux II 
: P /1 (m*p) _1 P 1/ : - : a 



[IV ' 



■ n 



//// 



and the stress tensor a 



(j,p 



The viscous part of the stress tensor, a' n = — a' 22 = 
-jPi^xVy + ^yVx), o' 12 = <r' 21 = \p{V x v x - V y v y ) is the 
Lorentz shear force. It is traceless, hence conservative. 



The effect can be interpreted in terms of semiclassical 
motion of electrons. A motion of electrons consists of a 
fast motion along small orbits and a slow motion of or- 
bits. A shear flow strains orbits elongating them normal 
to the shear, boundaries and vortices. Elongation yields 
an addition to the Lorentz, the Lorentz shear force which 
acts normal to the shear and proportional to the shear. 

An important consequence of the Lorentz shear force 
is accumulation of charges on boundaries and vortices - 
overshoots. They govern dynamics of edge modes [f| . 

12. Hall current. Another consequence of the Lorentz 
shear force is the increase of the Hall current by a non- 
uniform e.m. fields [2^. It follows from the conser- 
vation law (|25[) that in the linear approximation the 
Lorentz force equilibrates by the Lorentz shear force 
eE-iBv* w | Aw* fa |cA(J5/B). As the result 
the Hall current increases. In a uniform magnetic field 

3e = 0- xy (l + -t;(kl) 2 )E* k . 

So far we neglected diamagnetic effects by counting the 
energy from the ground state. The ground state energy 
depends on the magnetic field as —hMHu) c p, where M is 
an orbital moment per particle. In a non-uniform den- 
sity it yields the diamagnetic current jdia — — — MV*p. 
The diamagnetic currents gives a similar contribution to 
the Hall conductance as the Lorentz shear force. In- 
deed, the density of the chiral flow is determined by 
the flux according to the Eq. (fT7| . In the leading or- 
der m*Av ss 2nq'V*p ps mA(cE*/B), hence jdia ~ 
—o xy ^Mt 2 AE*. Thus diamagnetic current merely 

changes 4 > 4- + —M in the formula for the current. 

This result has been obtained in [22j where m*was set 
to be equal to bare electronic mass m e . A factor m»/m e 
offers an avenue to obtain the inertia m* by measuring of 
the increase of the Hall current by a non-uniform electric 
or magnetic fields. Comparing it with an independently 
measured gap allows to check Feynman-Bijl formula i 
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